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ABSTRACT 


A methodology for solar energy system design improvement has been 
developed and coupled with a constrained function optimization code 
resulting in an automated solar energy system design procedure. The 
scope of the methodology is limited to systems using flat plate 
collectors and water as the working fluid. 

Eight parameters have been included as independent design variables. 
The design variables included collector area, collector tile angle, 
collector and storage fluid stream velocities, and collector to storage 
heat exchanger dimensions. The procedure includes an accounting for 
economic parameters as an intimate part of the design process. The 
resulting methodology has been used for the design of solar energy 
systems which would use shelf item collectors for the purposes of 
determining the optimm design variable vector for a given situation. 
The methodology could also be used on a limited basis for collector design 
optimization by exploring the effects of changing selected collector 
parameters on system performance. The methodology is coded in the 
FORTRAN computer language under the name SOLOAD-1 (SOLAR ENERGY 
OPTIMIZATION ANALYSIS OR DESIGN). 

Initial system trials indicate complete stability with minimal 
constraint activations. Based on the results of approximately fifty 
design experiments using SOLOAD-1, new findings concerning optimum 
collector tilt angle and an invariant optimm collector flow factor 
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I. INTRODUCTION 


A. BACKGROUND 

The control and utilization of energy has become the 
major issue of this decade: the energy crises. The infla- 
tion rate of energy costs and particularly fossil fuel costs 
has the spector of uncontrollability as this decade comes to 
a close. Reference [1] reports a 30% increase in the price 
of oil in 1978 alone. Reports from the media in early Septem- 
ber 1979 indicate the price of home heating oil has incresed 
70% since the 1978 heating season. The crisis center appears 
to be an emerging realization that the supply of available 
energy is in fact exhaustible. This realization has spread 
from a few to the masses; unquestionably the fruit of the 
Organization of Petroleum Exporting Countries. 

In June of 1979, President Carter proposed to the Congress 
an energy strategy whereby solar energy would be providing 
20% of the Nation's energy by the year 2000. For the pur- 
poses of domestic hot water (DHW) and résidential space 
heating this goal is technically achievable. The basic solar 
energy technical theory is well documented by Duffie and 
Beckman [2] and Kreith and Kreider [3] and the technology 
is continually being updated in the Solar Energy Journal [4]. 
The increased effort in solar energy research work is in 
clear evidence as the documentation in this journal appears 


to be on an exponential rise. 
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The treatment of solar energy system performance is 
typically a determination of the fraction of a given heat 
load which is provided by solar energy. This fraction, f, 
is a nonlinear function of many system variables but typically 
increases in a monotonic fashion as a function of collector 
area. To simply seek maximum f and maximum collection of 
energy in an unconstrained fashion can result in extreme 
initial system costs. Accordingly, in order to assure a 
proper balance, technical performance and cost should be 
considered throughout the design process and economics should 
become an intimate part of the technical design problem; 
not an after thought at the end of system design. 

The economic ingredient in the problem of solar energy 
system design has received little attention in the literature. 
In view of this there is little information in the literature 
on how to proceed with design improvement of a solar energy 
system. Accordingly, a broad objective for this work was 
to develop a solar energy system model including economic 
considerations and to seek system design improvement by 


using the resulting eco-technical model together with an 


optimization algorithm. 


B. SCOPE 

The conversion of solar energy into useful work covers 
a broad spectrum of collection schemes. These include thermal, 
photovoltaic, biological, wind, and ocean thermal energy 


conversion. Reference [5] is an excellent introductory 


L3 


information source covering many methods of solar energy 
conversion and has become a classic. This thesis effort 

was limited to the thermal conversion of solar radiation 

for the purpose of generating low temperature heat. Since 
low temperatures are involved, only flat plate collectors 
were studied. Flat plat collectors can provide temperatures 
on the order of 200° F compared to concentrating collectors 
which can provide temperatures as high as 6700° F. The 
restrictions on temperature were chosen in keeping with 
availability analysis; matching the solar conversion system 
to its task. The task of providing high entropy, low quality 
DHW and space heating energy is matched with low availability 
(i1.e€., Minimum utilization of available energy in high 
quality, low entropy forms such as oil, gas, and central 
power). The collector working fluid was restricted to water. 
Further, since the collection of solar energy is a stochas- 
tic process, only analysis for long term performance was 
Studied. The analysis of the dynamic performance of specific 
systems in response to hourly climatic data has previously 
been conducted by Kline et al [6] and now forms the basis 


for long term analysis. 


fe OBJECTIVE 

The objective of this thesis was to develop methodolo- 
gies for: 

1. System synthesis including a determination of optimum 


system design variables including collector area, collector 


14 


' anniten Yiten palsevon 

qT .okeeniw = gmooed san boa a 

os 50 AOisvevrect Jemend ai of 
Lucsteqme?s wol gmaigeienep to 

5 io +212 yiae ,Bevioveal ous sete 

me shiverg nao smeteetion defq call bee 
polite 7a [80D Of Deteqeep F *o05 to 

I .4 “O05 fiplid oe sauccaisgmes sh ivoug a 
tw vo% ot noeonn ote oTueszegMes no SAG 


1o;e2evnon wiee_-«ts goisedam sieviaas) ve 


fs ea] seastime Aedes poebivony Zo Mealy ecT « 
dviw bedeted ef vptene paisset sosg 

ch yesees sldsliavs te nobtasstion erat oe 

isto bos ,ese ,llo ee dove wees ydoasne. wok: 


ot becSit>as: ssw biel? entdzow yostosLles eff 


jr aArrO8 iso oftes. viel 26? eleyvians vino « 
=u [¢ St remrtolted nimeryb ede Jo eleyians eff) 


Ss sivtamitio ylasnd 03 sanogee: £ 


ed oes 27 wet bite tar de @e acil”A ve becouks 


| 


on -aleylsne mxie3 9 


gviITt 
~loborgem golevab of Sew aeeey aids ‘So ‘svisoote 


umtge to seiterkinesel = pabhelom, wieadinya met ve 
ia x 


tarcallon ,se1s tosdelfos pasion eeidaivey mp tee: 
. wes a 2 = a Sar s+. 


wy 
ee a ee a ee F 


tilt angle, heat exchanger sizes, and flow rates. This 
analysis would be based on using shelf item collectors 
characterized by performance parameters Ft and FUL 
which have been determined by tests conducted in accordance 
with refs. [7 and 8]. 

2. Collector design optimization including a determin- 
ation of the optimum combination of collector design 
variables such as materials and geometries; and which 
would yield the optimum system performance. 

Fundamental to the effort was the utilization of: 

= COPES/CONMIN (Control Program for Engineering 
Synthesis with Constrained Function Minimization), a design 
improvement algorithm developed by Vanderplaats [9, 10, ll 
amd 12) . 

2. Standardized long term solar energy load fraction 
correslations (f - CHART) developed by Kline [6]. 

3. Economic considerations combined with technical 
analysis. 

4. Long term climatic data obtained from (NOAA) 


National Oceanic and Atmospheric Administration [13]. 
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IIT. SYSTEM MODEL 


A. OVERALL SYSTEM DESCRIPTION 

A schematic diagram of the solar heating system is 
shown in Figure 1. The system consists of four heat 
exchangers, a storage tank and associated system piping. 
The first heat exchanger is called the collector and is 
used to transform incident solar radiation into useful 
thermal energy. At the second heat exchanger the collected 
energy is transferred to the secondary loop and stored in 
the form of sensible energy in the storage tank. The 
sensible energy is transferred from storage to the loads 
via the third and fourth heat exchangers. Energy transfer 
is achieved via four forced convection flow loops. This 
model was further simplified by assuming total utilization 
of the energy stored in the tank and thereby eliminating 


the need to detail the DHW and space heating heat exchangers. 


B. ENERGY BALANCE 


The steady state power balance for the system is 
ore oF Wo aah 8, = eo (1) 


where, 


Useful rate of energy collection (Btu/hr) 


2) 
II 


W = Work energy utilized to maintain system flow 
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Q = Rate of utilization of auxiliary energy 


zs which is required in addition to solar 
energy in order to meet the DHW and space 
heating load 
Q, = Rate of energy loss by the system 


For the purposes of this thesis, it has been assumed 
that We is a known or estimated parameter. Most analyses 
in the literature do not account for W,, sance lt is typically 
assumed to be a very small fraction of the energy required. 
The percentage of the energy load which is supplied 


by solar energy is defined as: 


where the summation is taken over 12 monthly values. 

£ is a monotonically increasing function of collector 
area. No optimum f is evident. Figure 2 depicts a typical 
f distribution resulting from a design problem. Since 
initial costs are strongly dependent on collector area (see 
Section III), attempts to seek maximum f can result in very 
Midis initial system costs.. At.this point the motivation. to 


include economic considerations as a part of the system 


design model become clear. 


C. ECONOMIC ANALYSIS 
The matter of economics has been treated as it evolved 
in this study: a simple common sense approach. This 


approach includes: 
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1. An economic equation which utilizes technical 
parameters obtained from mass, momentum and energy balance 
considerations together with economic parameters. 

2. Standardization of the economic equation by char- 
acterization of all terms in present dollar values. 

The economic analysis is similar to an energy balance. 
Considering the economic system as the monetary account 
of a consumer and the time period as the economic lifetime 


of the energy conversion equipment: 


= - 
Din Pout D storage ie) 
where, 

Din = The value if dollars of the energy 
produced by the equipment 

D = The expense dollars to obtain, operate 

out ; ; : 
and maintain the equipment 

D = The net savings 

storage 


The rational consumer seeks maximum D This 


storage’ 
methodology is a simple capital budgeting technique and is 


contained in any standard reference in accounting or finance. 


oS “THE OBJECTIVE. FUNCTION 

Implicit in the use of the economic model is an assump- 
tion that the terms in the equation can be expressed as 
variables of the technical functions which have been 


defined by the physics of the engineering problem. For 
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the solar energy system utilized for domestic hot water 


and space heating the transformation is simple: 


Net Savings = Fuel Savings - System Costs 


This equation becomes the objective function by 


standardizing all dollar amounts to present values. The 


resulting equation is: 


= £ to 
NPV fQ.C Fi A.C, (4) 
where, 
NPV = Net present value of the solar investment, $ 
Ce = Fuel cost for the energy replaced by solar 


energy, $/Btu 


Bi = Present worth factor which standardizes fuel 
savings during system life into present 
dollars, years 


A = Collector area, eee 


€ = System lifecycle cost per unit collector 
area, $7 £2 
The effort now proceeds to determining each of the 
above dependent parameters and to obtaining those combina- 
tions of these variables which will result in maximum 


values of NPV. 
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III. SYSTEM OF EQUATIONS 


A. GENERAL 

The objective function has been fully developed by 
examining each of its parameters individually. Completion 
of this development has resulted in a non linear function 
of eight independent design variables; a hypersurface in 
eight space. When formulated in this manner, the problem 
is clearly not amenable to analytical solution. Ina very 
simple manner however the problem can be interfaced with 
the design improvement algorithm COPES/CONMIN (see Section IV). 
The primary contribution for this effort comes from Kline [6] 
who developed the f-chart correlations for determining monthly 
solar energy load fractions fi The annual load fraction 
£ is then computed by a weighted average of the monthly 


values. 


Bs. SPECIFIC PARAMETERS 
1. Annual Energy ,Load, QF 

The annual energy load consists of two basic ingredi- 
ents; the DHW load and the space heating load. Assuming an 
average mean ground temperature the month to month DHW load 
was considered constant except for the variation due to month 
length. The methods for determining the space heating loads 
typically follow the ASHRAE manual, reference [14]. Refer- 


ence [15] is recommended as a text for space heating load 
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computations. For the purposes of this study representative 
building loss parameters were used. Typical building con- 
ductance or (UA). values of 30,000, 20,000 and 10,000 
(Btu/deg F day) have been used based on a standard building 
of 1750 £t? of floor area, a heat transfer surface of 5000 
re and building loss coefficients, ULs OE 625 elle, Ane «09 
Btu /ft- hr F respectively. The conductance, UA, is the 
space heating load at design conditions which has been 
estimated in the manner of Ref. [14] and divided by the 
design temperature difference. 
2. Annual Solar Energy Load Fraction, f 

f is a function of collector performance parameters, 
thermal physical properties of the working fluids, fluid 
flow rates, heat exchanger performance parameters, collector 
falc angle, climatic conditions and latitude. The climatic 
conditions determine the load distribution (heating degree 
days), the solar energy flux distribution (insolation) and 
the environmental stress on the collector (ambient tempera- 
Eure). 

a. Collector Performance 

A determination of collector performance proceeds 

from the well known collector equation of Hottel and Whillier 


[16] which is the result of an energy balance on the collector: 


om = PUULAL(IRta - U, (tT; - eee (5) 


zl. 


os 
i 


where, 


= = Useful rate of energy collection (Btu/hr) 

e = Collector heat removal factor (dimensionless) 

UL = Collector loss coefficient (Btu/hr £t* F) 

A, = Collector area as 

f= solar insolation on a horizontal surface (Btu/hr Ft7) 

R = Collector tilt correction factor (dimensionless) 

T = Transmittance of collector cover system (dimensionless) 

a = Absorptivity of the collector absorber plate 
(dimensionless) 

T. = Collector fluid inlet temperature (F) 


Ta = Climatic ambient temperature (F) 

The analytical expressions for ee and Uy, are complex 
and result from a lengthy development which can be found in 
references [2 and 3]. The computation of Ps and UL using 
these expressions is not necessary for system design analysis 
using shelf item collectors since the parameters FUL and 


Fto can be obtained from collector efficiency tests. 


Collector efficiency is defined as: 


7 = (6) 


Assuming that U, is not sensitive to changes 


L 


in the environment and combining equations (4) and (5) results 


ans 


Equation (7) is linear in the collector 
parameter (T, - T,)7IR and forms the basis for determining 
collector performance. It is clear that the parameter 
FUL is obtained from the slope and FTO is obtained from 
the intercept of equation (6). Some typical curves of 
collector efficiency are included in Figure 3. 

Two standards have been developed for collector 
testing to measure collector performance, namely references 
f71 and [8]. The standards differ in the independent 
variables used in the performance characterization. 
Reference [7] the ASHRAE standard uses T = T. = Toi 
ref [8] the NBS standard uses T = (T. + T,) 72. where t, 
is the collector outlet temperature. All collector 
performance test data used in this effort is based on the 
ASHRAE standard. All collector data used in this effort 
is based on actual collector tests conducted by NAVFAC 
and reported via ref [17]. 

Be Collector to Storage Coupling 

The collector becomes coupled to’the storage 
tank via the collector to storage heat exchanger following 
the method of de Winter [20]. It is assumed that the rate 


of energy transfer at the heat exchanger is equal to the 


rate of useful energy collection: 
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T = Maximum system temperature or collector outlet 


Hy 
iT} 


Minimum system temperature or storage outlet. 


Equations (5) and (8) are combined to eliminate 


T; and ed dependence to give: 


a = PIU, (IRta ~ Sed - Eee, (9) 
where, 

E = Bex (10) 

ee (mc_) = 
oe = {1 fe See ae Pe 11} (11) 

ex . . 
(mc_) etme)... 
p p min 


eonncolLagyinsolation,at Optimum TiLt Angle, Ty 
The solar radiation intensity on.a tilted 


surface is by definition: 


Solar insolation is typically measured and 
reported for a horizontal surface. Improvements in a 
sOlar collector installation are sought by tilting the 


collector to the optimum angle. The computation for 
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optimum angle is iterative for each design situation. 
Several rules of thumb have evolved: 

(1) Collectors should be oriented at a slope 
of .9 times the latitude angle for maximum annual collection. 

(2) For DHW systems where loads vary little 
during the year the best angle of tilt is equal to the 
Petude angle. 

(3), (fhe: optimum collector tilt angle for least 
cost per Btu delivered for building heating is approximately 
the latitude angle plus 15 degrees. 

The methodology for determining collector tilt 
angle used in this effort was as follows: 

(1) Horizontal monthly data was obtained for 
97 locations from NOAA in reference [13]. 

(2) An algorithm was developed to compute 
monthly horizontal extraterrestrial radiation intensities. 
fais computation is: a function of latitude, daily hour 
angles on horizontal and tilted surfaces, daily declination 
angle, and collector tilt angle. The azimuth angle was 
always chosen for a due south collector orientation. 

(3) The methodology follows that of Liu and 


Jordan in reference [19]: 


R f=*- (r= 
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Go = 1.3903 - 4.0273(4) + 5.5315(4)7 - 3.108(2)7 (14) 
fe) fe) fe) 

Ww. ‘Sin(o—s)sin_6_+ cos({o—-s)cos 6 sin w_' 

Re = eo ess: (15) 

d w sin iesim oO + Coss @ cos 6 sin W, 

. : ale sil) = 
oS = min[cos ~(-tan(o-s)tan 6), cos “(-tan » tan 6)] (16) 
=i 

ws = cos [-tan ¢ tan 6] (17) 

R aoe (i +.cos s) (18) 

d 2 
1 

> (l - cos s)op (19) 

where, 

R = monthly slope correction factor 

Rp = average ratio of monthly direct radiation on a 
Eiited surtace to monthly direct radiation on a 
horizontal surface 

R= ratio of monthly diffuse radiation on a tilted 
surface to monthly diffuse radiation on a 
horizontal surface 

R = xratio of monthly reflected radiation on a tilted 

P surface to monthly total radiation on a horizontal 

surface 

¢ =. ratio.of monthly diffuse radiation on a horizontal 
surface to monthly total radiation on a horizontal 
surface 

- = ratio of monthly total radiation on a horizontal 

fe) surface to monthly total extraterrestrial radiation 

On a horizontal surface 

LF Sunrise hour angle on a horizontal surface 

p= Sunrise hour angle on a tilted surface 
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s = collector tilt angle 


¢ = latitude angle at collector site 
6 = declination angle of the earth 
9 = ground reflectance for area adjacent to collector. 


3. Fuel Cost Parameters 
Fuel cost parameters included a selection from among 
meruiet! types (i.e. oil, electricity, or gas), unit of issue 
cost, fuel heating value and efficiency of the auxiliary 
heating system. The SOLOAD-1 system allows complete freedom 
in the selection of these parameters. A typical set of fuel 


cost parameters used in a design experiment included: 


Fuel: Oil 

Unit Cost: 9 ($/Gal) 

Heating value: 142,000 (Btu/Gal) 
ELELCLency: aif! 

Resulting Fuel Cost: Ce = $9.05 per Tokg Btu 


4. System Cost Parameters 


System cost considerations typically include 
initial costs and life cycle costs for operation and 
Maintenance. Detailed guidance for complete cost consi- 
derations is included in ref. [3]. Further, reference [20] 
contains guidance to be used in determining costs for US 
Navy installations. The following cost parameters were 


used in SOLOAD-1: 
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where, 
CL = Total system life cycle cost 
- = Total cost per square foot 
Cc = Cea (27°) 
Ss a Vv 
= - - ~ 
Cc; (sy CL Che c.. (22) 
S = OP (C rs FC.) (2.3,) 
Vv ai ea 
where, 
C; = Initial system cost 
Cc. = Operation and maintenance costs 
Ce = Collector cost 
Cr = System installation cost 
a Storage tank costs 
C.. = Collector to storage heat exchanger costs 
co = Pumping power costs 
ei = Maintenance cost factor 
Fi = Present worth factor 
yy oN 
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where, 


1 eerie: 
a (25) 
i* = Discount rate 
j = Fuel inflation rate 
N = System life in years 


It should be noted that Fi occurs explicitly in 
the first term of the objective function and implicitly 
in the second term. It should also be noted that initial 


costs are not amortized in the SOLOAD-1 model. 
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IV. NUMERICAL OPTIMIZATION 


A. SIMPLE OPTIMIZATION 

Design problems typically seek the minimization or 
maximization of an appropriate parameter within a frame- 
work of constraint specifications. The parameter to be 
optimized may be a function of several design variables and 
is termed the objective function. Other parameters which 
may be separate functions of design variables must not 
exceed specified bounds for the design to be acceptable. 
These parameters are termed design constraints and are not 
to be confused with limits which may be set on design 
variables which are usually termed side constraints. 

Engineering problems can be numerically coded for an 
analysis (once through) solution. The simplest scheme for 
optimization may consist of a series of loops through the 
computer code which may cycle many combinations of design 
variables. The combination of variables which provides the 
best design and which also satisfies the constraints is 
then considered optimum. This approach may be economical 
for a design problem with just a few design variables and 
short computer time requirements. A design problem with 3 
design variables, ten values for each design variable, and 
One-tenth seconds central processing unit (CPU) time for 
each analysis would take a total of 100 seconds of CPU time. 


The solar energy optimization design problem as characterized 
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by the code developed by this thesis has a minimum of 8 
design variables; each analysis of its objective function 
requires about 2 CPU seconds. Using the simple approach 
and assuming ten values for each design variable would 
result in a CPU time of 68 years for each design problem. 
Clearly, a more rational approach to optimization is 
necessary. 

Vanderplaats [9] suggested that many special algorithms 
for numerical optimization have been proposed in recent 
years, but that in many cases unsuspecting practitioners 
find their particular optimization problem unsolved only 
after large amounts of time and effort are expended. This 
can occur usually because of inexperience by the practitioner 
in determining the limitations of specified algorithms. 
Vanderplaats [10, 11, and 12] has developed a FORTRAN 
coded algorithm capable of optimizing a very wide class of 
Engineering problems. The system includes COPES (Control 
Program for Engineering Synthesis) and CONMIN (Constrained 
eunceron Minimization). “This’ optimization ‘System is referred 


to as COPES/CONMIN. 


B. THE COPES/CONMIN SYSTEM 


1. Terminology 
CONMIN is a FORTRAN program in subroutine form for 


the solution of linear or non-linear constrainted optimiza- 


tion problems. The user prepares an analysis program. The 


oui 


program must be named SUBROUTINE ANALIZ. The process of 
computer aided design or of trade off studies with a minimum 
of man-machine interaction becomes fully automated via the 
COPES program. Three basic definitions are required: 

Design Variables. Those parameters which the 
optimization program (CONMIN) is allowed to change in order 
to improve the design. Design variables appear only on the 
right hand side of equations in the analysis program (ANALIZ). 
Limits imposed on design variables are termed side constraints. 

Objective Function. Usually the single parameter 
which is to be minimized or maximized during optimization. 
The objective function always occurs on the left side of 
the equation in the analysis program. (Refer to [12] for 
exceptions.) The equation defining the objective function 
may be linear or non-linear, implicit or explicit, but must 
be a function of the design variables to be meaningful. 

Design Constraints. Any parameter which must not 
exceed specified bounds for the design to be acceptable. 
Constraint parameters always appear on the left side of the 
constraint function equations. Constraint functions may be 
finear or non-linear, implicit or explicit, but must be 
functions of the design variables. 

Assuming that the optimization process requires the 
maximization of a particular objective function, the general 


Optimization problem can be stated as: 
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a. Find: X which maximizes f (X) 
be “sub Jece: to: 
(1) CONSTRAINT EQUATIONS G(X) Uy, 3) = abe NEON 


(2) SIDE CONSTRAINTS VLB cS x. < VUB., 1 = 1, NDV 


Where, X = X(X)+ X5, +--+, X,) is the vector 
of NDV design variables, F(X) is the objective function and 
G(X) are the set of NCON constraints. VLB, and VUB, are 
the lower and upper bounds respectively on each design 
variable. 

2. Methodology 

The solution process proceeds as follows: 

a. The user prepares an analysis subroutine which 
defines X, F(X) and G(X). This subroutine must be named 
ANALIZ. ANALIZ must have three segments; input, analysis 
ano Output keyed to COPES flags, ICALC = l, 2, or 3 
respectively. 

b. The user prepares an input data file for COPES 
which includes a wide variety of system options, appro- 
priate matching mechanisms between ANALIZ and CONMIN and 
the constraint boundaries. 

c. Using the initial vector of design variables 
COPES obtains an initial solution from ANALIZ and subsequent 
solutions by updating X as determined by CONMIN. Any analy- 
Sis sOlution which satisfies the constraint equations and 
the side constraints is a feasible design. If an analyti- 
cal solution violates any of these constraints it is an 


infeasible design. The minimum feasible design is optimal. 
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The feasibility determination includes: 

(iy Liva constraint equation is violated 
ooo. if G(X) > 0) then the jth constraint is violated. 

(2) If a constraint equation equality 
Gendition is met (i.e., if G(X) = 0) then the jth con- 
straint is active. 

(3) If a constraint equation condition 
ms met (i.e., if G(X) < 0) then the constraint is inactive. 

Note that CONMIN is designed to minimize objec- 
tive functions; the process of maximizing an objective 
function is concerned with minimizing the negative of an 
objective function. 

(4) Similarly, side constraints may be 
inactive or active but side constraints will never be 
violated in a particular analysis computation because 
they are specified limits not dependent variables which 
Baetne: case for design constraints. 

(5) All inequality conditions are repre- 
sented by a band around the zero condition due to computer 
limitations in defining zero. 

Ge ot the 1nitral analysis: solution is? infeasible 
CONMIN moves towards a feasible solution by adjusting the 
design variables appropriately. The optimization process 
then proceeds in an iterative fashion as follows: 

(1) The iterations are governed by the 


recursion relation: 
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where, 
q = iteration number 
a* = a scalar move parameter which defines the 
distance of travel in the direction of 
search 
S$ = direction of search 
(2) S is determined such that F(X) will 
be minimized without violating any constraints. CONMIN 


calculates the gradient of F(X) by using finite difference 
techniques. Because no constraints are violated, the 
greatest improvement in F(X) will be realized by moving 

in the direction of steepest ascent, the gradient of F(X). 

(3) Once the directive is known, the 
move parameter, a*, which will allow the largest magnitude 
improvement in F(X) is to be found. A one dimensioned 
search of the S direction is carried out until the best 
improvement number, a*, is found. 

(4) CONMIN utilizes methods other than 
the method of steepest descent in determining S particularly 
in the presence of active constraints. These methods 
include the method of conjugate directions developed by 
Fletcher and Reeves [21] and the method of feasible direc- 


tions developed by Zoutendijh [22] and implemented by 


Vanderplaats and Moses [23]. 
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e. CONMIN continues to iterate for an optimal 
design by computing successive S and a* values always 
keeping within the defined constraints. If there is no 
relative or absolute change in F(X) for three successive 


iterations, the optimum is considered found. 
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V. “SOLOAD=1 SYSTEM 


A FORTRAN coded algorithm for the analysis of systems 
containing shelf item collectors has been developed. The 
algorithm has been named SOLAR ENERGY OPTIMIZATION ANALYSIS 
OR DESIGN (SOLOAD-1). The portions of the algorithm 
reserved for design optimization of collectors were not 
completed and have not been included in SOLOAD-1l. The 
interface between SOLOAD-1 and COPES/CONMIN has been 
developed via the Control Program (CP-67) and the Cambridge 
Monitor System (CMS) a time-sharing system developed for 
the IBM System 360 Model 67. SOLOAD-1 includes: 

A. Three executive programs to provide the interfaces 
among the system elements and to initiate design problems. 

B. An analysis subroutine which includes the system 
of equations covered in Section III. A summary of the 
objective function, design variables, and design constraints 
is included in Appendix A. 

C. Three auxiliary subroutines including a calendar 
data array, climatic data and user defined input parameters. 
D. Data files for user input to SOLOAD-1 and COPES/ 

CONMIN. 

The interface among SOLOAD-1 system elements and 
COPES/CONMIN is shown in Figure 3. SOLOAD-l programs and 
Subroutines are included in Appendix B. Data files are 


included in Appendix C. 


FP 


Vii. « RESULTS 


A. GENERAL 

The first objective of developing a methodology for 
the analysis of systems with previously defined collector 
design was achieved. The second objective of developing 
a methodology for optimizing the collector component design 
was not achieved. However, it appears that the potential 
exists for using the system analysis methodology to obtain 
improvements in collector design. For example, if cost 
data on various types of collector cover plates or collector 
absorber surfaces can be correlated with collector per- 
formance parameters, then the system model could be used 


for quick checks on collector design improvements. 


B. VERIFICATION 

A single design problem was selected and started using 
10 different sets of the starting design variable vector. 
The same optimal design result within + 0.5% was achieved 
for each run. This was the only actual verification deemed 
possible at this time due to lack of any known optimization 


data in the literature. 


C. EXPERIMENTS 
Approximately 50 design experiments were analyzed. 
Each problem was characterized by a unique identification 


number which could tie together the location, economic 
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environment, the collector, and the space heat transfer 
coefficient, UA, for the experiment. For each experiment 
an input parameter summary report and output parameter and 
results summary report was generated. These reports are 


included in Appendix D. 


D. SPECIFIC RESULTS 
te Constraints 
The only active design constraint in most experi- 
ments was the tube thickness in the collector loop at the 
heat exchanger. Some experiments resulted in no active 
design constraints. There were no active side constraints 
in any experiment. 
Zea Collector, Flow: Factor 
The collector flow factor F" is defined as the 
ratio of the collector heat removal factor to the collector 


Serr eLvency factor: 


ny 


’ 1g 
ay = rae 


The flow factor result obtained in all optimization experi- 
ments was the same. F" = .948 + .008. 
ee Collector Tilt Angle 
Each experiment conducted for a particular geo- 
graphical location resulted in the same collector tilt 
angle. However, experiments for different locations but 
the same latitude resulted in significantly different 


optimum collector tilt angles. 
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VII. CONCLUSIONS 


A. The SOLOAD-1 system in conjunction with COPES/CONMIN 
appears to offer the potential for further improvement and 
potentially a valuable automated technique for solar energy 
system design. 

eeerke “rule of thumb"=stypically used for»collector tilt 
angle optimization (i.e., latitude plus 15°) should be 
used with caution since preliminary results indicate a 
strong tendency for climatic dependency. 

C. The uniqueness of the flow factor was suggested by 

the continued result of .948 + .008 for each experiment. 
Pending confirmation by further testing it appears that 

a simple correlation for determining optimum collector 
loop flow rates may be available. Based on a flow factor 
Of 0.948 the resulting correlation would be: 


Gypse 601955 PU (gpm/£t* of collector area) 


L 


This correlation is developed in Appendix E. 
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VIII. RECOMMENDATIONS 


A. The basic model of SOLOAD-1 should be upgraded to 
include specific characterization instead of simple 
parameter selection and input for: 

1. Pumping power in all four loops 

2. Inclusion of the DHW heat exchanger 

3. Inclusion of the space heat exchanger 

4. Inclusion of building loss coefficients as a 

design variable. 

B. Prior to any additional experiments with the present 
model, a complete survey of the industry should be con- 
ducted for collector parameters and costs. Following this 
survey, a series of experiments should be conducted in 
search of correaltions among collector performance parameters 
Pie. ; FL (ta) and FUL), collector unit costs, and system 
performance. 
C. A larger sample size should be used to verify the 
uniqueness of the optimum collector flow factor (F") as 
Suggested by the results of 50 experiments. 
D. The remainder of the NOAA climatic data bank (i.e., 67 
more cities) should be included in SOLOAD-1l. The optimum 
collector tilt angle could then: be computed for the 97 
NOAA regions. 
E. The model should be upgraded to accommodate analysis 


for systems with air as the working fluid. 
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F. The model should be upgraded to accommodate amorti- 
Saetion Of the installation cost instead of just initial 


cash payment. 
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APPENDIX A 


SUBROUTINE ANALIZ SUMMARY 
A. OBJECTIVE FUNCTION 


= = Y= 
NPV £ Qr Ce Ps ie or 


B. DESIGN VARIABLES 

| Collector Area A 
Collector tile WAngle cS 
ColLtlector,Loop«I.D. d 
Collector Loop’ 0.D. d 
Heat Exchanger I.D. d 
Collector Flow Velocity ve 
Storage Flow Velocity Vv 


Heat Exchanger Length L 


C. DESIGN CONSTRAINTS 


G) = qe = dq. 

oP = dg - o 

G, = Reynolds number, Collector loop 

G, = Reynolds number, Storage loop 

Ge =" Capacity ratio, A om 

Ge = Ist Flow parameter, a = Be ee Us, 
Ga = 2nd Flow parameter, c> = Be ee a 
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APPENDIX B 


SOLOAD-1 Computer Program 
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APPENDIX D 


EXPERIMENT REPORT SUMMARIES 


Each optimization problem or experiment is represented 
by an input summary report and an output summary report. 
These reports are tied together by a unique identification 


as follows: 


where, 


LL = location identification number in 
accordance with Appendix C 


Ny = present worth factor identification number 
ree = LB aoe 
cu 
Z2=Pr = 24.34 
i 
N, = collector identification number 


1 - Solarnetics 
2 - American Sun 


3 - Federal Prison Institute Double Glaze 


N, = heat load conductance 
1 - UA = 30000 Btu/HDD 
2 - UA = 20000 Btu/HDD 
3 - UA = 10000 Btu/HDD 
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APPENDIX E 


POTENTIAL CORRELATION FOR OPTIMUM COLLECTOR FLOW RATE 


Reference [2] defines the collector heat removal factor, 


Ge -F'u, 
a = Te {l - Gace, (Bl) 


Further, the collector flow factor, EF", is: defined as, 


F Ge “F'UL 
F = Sees ru i = exp ( Ge ) } (B2) 


Reference [2] shows that F" approaches the unity value 
asymptotically as the parameter Co eae increases. 
Reference [3] recommends capacity rates of 10-15 lb/hr =e 
(e0223*—= 0334 Gpm/ft®) as the best compromise among 
collector heat transfer coefficient, fluid pressure drop, 
and energy delivery. The unit ft. refers to collector area. 

The results of this thesis suggested a correlation 
between the collector performance parameter FUL and the 
Capacity rate. Once these results are verified by further 
testing including model changes to include fluid pressure 


drop, a simple correlation may be available and follows 


from: 
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Let, K* = The optimum flow factor determined from 
computer experiments. 


Ss es bone Uy 
oe 


fie follows from equations (Bl) and (B2) that: 


‘aad = = 
2 
zy = 1- exp (-1/z,) 
Or, 
ee ee ae. 
* 

Let, [en The optimum parameter which corresponds to 
K* and is obtained by solving equation (B2) 
for & 

aL 
Therefore, 
oe 
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Based on the results of a limited number of computer 
experiments all yielding an apparent invarient flow factor, 


F" = .948: 
eo = 201955 (Com hr F/BEu) 
Om, 


Z 
G = -01955F_U (Gpm/ft_) 
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APPENDIX E 



POTENTIAL CORRELATION FOR OPTIMUM COLLECTOR FLOW RATE 



Reference [2] defines the collector heat removal factor 
F r , as: 

Gc -F ' U 

F r = -o* {l-exp(-^-i)} (Bl) 

L p 

Further, the collector flow factor, F", is defined as, 

F Gc -F ' U 

F " = r- - f< (1 - -get 1 ” <b2) 

L p 

Reference [2] shows that F" approaches the unity value 

asymptotically as the parameter Gc /U F ' increases. 

P U 

2 

Reference [3] recommends capacity rates of 10-15 lb/hr ft 

2 

(.0223 - .0334 Gpm/ft^) as the best compromise among 
collector heat transfer coefficient, fluid pressure drop, 
and energy delivery. The unit ft refers to collector area 
The results of this thesis suggested a correlation 
between the collector performance parameter F^U^ and the 
capacity rate. Once these results are verified by further 
testing including model changes to include fluid pressure 
drop, a simple correlation may be available and follows 
from: 
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Let , 



k* = The optimum flow factor determined from 
computer experiments. 



=1 - Gc p /f ' u l 
?2 " Gc p /F r U L 



It follows from equations (Bl) and (B2) that: 



F" 



1 



5 



2 



■— = 1-expt-l/?,) 

^ 1 ± 



? 1 -ln(l - l/^ 2 ) 



Let, = The optimum parameter which corresponds to 

<* and is obtained by solving equation (B2) 
for c. 



Therefore , 



or 



or 



Gc /F U T 
p r L 



JL 

. ★ 
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! 



or. 



G = 



c i 



F U/c 
r L p 



or, 



G = k**F U t 
r L 



Based on the results of a limited number of computer 
experiments all yielding an apparent invarient flow factor, 
F" = .948: 



<** = .01955 (Gpm hr F/Btu) 



Or, 



G 



. 01955F U_ 
r L 



(Gpm/ f t^ ) 
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